Delay of Disorder by Diluted Polymers 
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We study the effect of diluted flexible polymers on a disordered capillary wave state. The waves 
are generated at an interface of a dyed water sugar solution and a low viscous silicon oil. This allows 
for a quantitative measurement of the spatio-temporal Fourier spectrum. The primary pattern after 
the first bifurcation from the flat interface are squares. With increasing driving strength we observe 
a melting of the square pattern. It is replaced by a weak turbulent cascade. The addition of a small 
amount of polymers to the water layer does not affect the critical acceleration but shifts the disorder 
transition to higher driving strenghs and the short wave length - high frequency fluctuations are 
suppressed. 



It is known from the 1950s 1] that small amounts of 
flexible polymers can cause drag reduction in turbulent 
flow. But both experimental and theoretical investiga- 
tions have not yet conclusively answered the questions 
what the fundamental mechanisms of drag reduction are 
and statistical and stability arguments coexist 2] . A dif- 
ficulty in the theoretical description of drag reduction 
in fully developed turbulence (e.g. in pipe flow) is that 
the polymer effect, the drag reduction, occurs in a range 
where any theoretical treatment of the flow must be fully 
nonlinear. So called weak turbulence |^, reported for dif- 
ferent systems like semiconductor lasers an waves, obeys 
similar properties like fully developed turbulence, e.g. 
the energy is injected on large scales and transported via 
a continuous cascade to shorter, viscous length scales. 
But typical Reynolds numbers are much lower and for 
the analytical treatment small background waves can be 
taken as the small parameter around which a linearisa- 
tion can be performed. 

In this Letter we present experimental measurements 
that show how polymers can affect the transition to a 
weakly turbulent state and delay the disordering. We 
present quantitative data of parametrically driven waves 
(Faraday waves, see Fig. CJ at the interface of two liquid 
layers. Due to the low viscosities the system evolves al- 
ready at low driving strengths from the square state into 
disorder. Our measurements show that by adding even 
of a small amount of flexible polymers (Polyethylenoxide, 
PEO) to the aqueous layer the transition into disordered 
state occurs at considerably higher driving accelerations 
whereas the critical acceleration of the primary bifurca- 
tion remains unaffected. We will see that the pronounced 
elongational viscosity of the polymer solution accounts 
for this behavior. 

The experimental setup is described in detail elsewhere 
In order to overcome the problems associated with 
light refraction at the interface of a free surface we chose 
to work with a system of two layers of refractive index 
matched unmiscible liquids of different densities. The re- 
placement of air as the upper fluid for free surface waves 




FIG. 1: Snapshots of the interfacial waves at a driving am- 
plitude ao = 31.7 m/s^. The lower fluid layer is (a) the pure 
dyed water sugar solution and (b) with 110 ppm PEO. (c) 
and (d) are the corresponding power spectra. 



is hydrodynamically nothing but an increase in viscosity 
and density, but the actual interfacial deformation can be 
simply determined by the absorption of light passing the 
lower, dyed layer. The method has already been proven 
to allow for quantitative determination of the complete 
spatio-temporal spectrum of ordered wave states, with a 
resolution in amplitude better than one percent of the 
maximum wave height. The lower liquid is a sugar- 
dye- water mixture (density p = 1200/c^/m^, viscosity 
T] = 3.2 mPas) in which the high molecular PEO poly- 
mer {Mw = 6 X lO^amu) can be dissolved. The highest 
polymer concentration is 110 ppm well below the critical 
overlap concentration c'^ = 300 ppm. The upper liquid 
is a silicon oil (density p = 910kg /m^^ viscosity r] = 0.65 
mPas). The low viscosity was chosen to minimize the 
influence of the upper layer, similar like in free surface 
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waves. The interfacial tension is cr = 35 dynes/cm, yield- 
ing a capillary length of Ic = (c^/ {^PoY^'^ = ^Tnm^ with g 
the earth acceleration and Ap the differences in the fluid 
densities. The thickness of the water sugar (silicon oil) 
layer is ho = 3mm (1mm). A closed circular vessel with 
glass plates at the top and bottom is attached to an elec- 
tromagnetic shaker. The temperature is hold constant at 
T = 23 ± 0.1°C. The system is accelerated in the form 
a{t) = aosin{Qt) and the driving frequency is always 
Q/27T = 47Hz. We are in the regime of the generic sub- 
harmonic response of the Faraday Experiment and the 
waves oscillate with a fundamental frequency coq = Q/2. 
It is related via the dispersion relation to a wavelength 
27r//co = 5mm and we are close to the infinite depth limit 
tanh{koh) = 1. Our measurements are performed by in- 
creasing the driving strength ao in steps of 2%. After an 
equilibration time of 60 seconds the interface is filmed for 
two seconds by an 8-bit camera with a spatial resolution 
of 256 X 256 pixels at repetition rate of 500 Hz. The max- 
imum acceleration is ao = 95m/s^ and for ao > 54m/5^ 
first droplet ejections occur, but only for ao > 70m/ s'^ 
the interface disintegrates significantly. 

All solutions show the same critical acceleration ac = 
17.5 m/s^ at which the interface becomes unstable first. 
Near onset a pattern of squares is formed. Figuredshows 
a typical snapshot of a measurement at a driving am- 
plitude of ao = 31.7 m/s^ with the lower layer (a) the 
water sugar solution and (b) after the addition of 110 
ppm PEO. Strikingly, for the diluted polymer solution 
the quadratic order is still clearly to identify while the 
order in the water sugar solution is already " melted" . A 
typical temporal realization of the interface displacement 
h{x^ t) is shown in Fig. [21 For the polymer solution we 
observe a small drift while the water sugar solution is 
strongly intermittent. 

The melting processus can be characterized by the nor- 
malized angular autocorrelation function C(0°), defined 
by van de Water and Binks 6] such that square pattern 
will lead to a pronounced peak at C(0 = 90°) (Fig. [2|). 
For the polymer solution a quadratic order is found up to 
driving strength ao ~ 30m/5^. For the water sugar solu- 
tion the C{(j) = 90°) peak is much weaker and disappears 
for ao > 20m/5^. 

Figure shows the temporal and angular averaged am- 
plitudes of the Fourier spectrum. The prominent 
peak at the basic wave number ko is clearly recognized. 
This is the length scale at which energy is injected into 
the system. However, already at intermediate driving 
strength ao = 31.7 m/s^ the higher harmonic spatial 
Fourier peaks are very weak and the spectrum almost 
consists of a continuum only. The theory for weak tur- 
bulent capillary waves 3] predicts a wave number density 
equivalent to 

Fa ~ k-''/\ (1) 
The exponent could be confirmed by experiments with a 
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FIG. 2: A typical temporal evolution of the interfacial defor- 
mation at an arbitrary point near the center of the container. 
Parameters are like in Fig. Q 

free surface ^] and our interfacial waves show the same 
dependency at higher driving strengths, too (see straight 
line in Fig. 

We find a 5% difference in the basic wave vector ko for 
the water and the polymer solution but for lower driving 
strengths the main difference appears at large wave num- 
bers. For the polymer solutions the fluctuations in the 
small length scales - equivalent to high temporal frequen- 
cies 7] - are almost one order of magnitude lower than for 
the solution without polymers. The integrated intensity 
F/ (see Fig. |5| is dominated by the power of the basic 
wave vector and is similar for both systems at ao < 40 
m/s^ F/. This is in contrast to earlier Faraday Exper- 
iments with more viscous polymer solutions j^, where 
e.g. the critical acceleration is significantly affected due 
to viscoelasticity. Only for driving accelerations ao > 40 
m/s^ the integrated intensities Fj of the solution with 
polymers is lower than for the pure solution (see Fig. 21 
and[5J. 

We think that the phenomena observed in our exper- 
iments can be only interpreted if one takes into account 
the differences in elongational and shear viscosity. A sim- 
ple, Newtonian liquid is defined by a constant viscos- 
ity T] that does not depend on the applied shear rates 
or flow forms. Already in linear response the viscos- 
ity of viscoelastic liquids becomes a complex quantity 
ri*{(jo) = Tj' {(jo) -\- ir]"{uj), with rj' the dissipation and rj" 
the elasticity. It can be measured as e.g. in an oscilla- 
tory shear experiment with small shear amplitudes. In 
a similar experiment with constant shear rate 7 polymer 
solutions might show a shear thinning behavior and 77(7) 
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FIG. 3: (Color online:) The temporal averaged and base- 
line subtracted angular autocorrelation functions (3(0°) of the 
Power spectrum of the interfacial deformations. The color 
density plot shows the evolution of the autocorrelation func- 
tion versus driving strength (left y-axes). The quadratic order 
in the water sugar solution is weak and a standard plot of the 
autocorrelation function at ao — 25m/s^ is shown as an ex- 
ample (black line, right y-axes). The lower fluid layer is (a) 
the pure dyed water sugar solution and (b) with 110 ppm 
PEG. 
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FIG. 4: (Golor online) The temporal and angular averaged 
squared amplitude Ya of the Fourier spectrum of the inter- 
facial states versus the radial averaged wave number k. Each 
curve is the average of 1000 images of the interfacial wave 
state. The lower liquid is the water sugar solution (symbols) 
or the 110 ppm PEG solution (lines). The driving amplitudes 
are ao = 31.7 m/s^ (triangles and red line) and ao = 60 m/s^ 
(squares and green dashed line). The short black line indi- 
cates a -15/4 slope. The arrow indicates the position of the 
first higher order peak at v^/co- 



depends on the shear rate. But shear thinning is more 
pronounced for concentrated solutions. The viscoelas- 
tic solutions used in the earlier Faraday Experiments |^ 
are characterized by a significant linear elastic contribu- 
tion T]" {(jj) and shear thinning. Their interaction with 
the Faraday waves and the influence of the linear stabil- 
ity could be well explained in terms of a linear response 
theory. 

The diluted low viscous polymer solution used in our 
experiments does not show an elastic contribution rj' {uo) 
that could be measured with a high resolution rheometer 
(MARS, Thermo Electron Germany, Karslruhe) and its 
shear viscosity measured at different constant shear rates 
1 < 7 < 100 differs only a few percent from the solvent 
viscosity. But diluted solutions of flexible polymers are 
known to possess a pronounced elongational viscosity r]e 
Any flow can be separated into a rotational and 
an elongational part. This refers to the symmetric and 
the antisymmetric parts of the velocity gradient tensor L, 
with components Lij = \/jVi {vi are the velocity compo- 
nents of the flow). For Newtonian liquids rje is just given 
by a simple geometrical factor, the Trouton ratio with 
r]e = 3r]. A pure rotation has no effect on the polymers 
but an elongation stretches them most efficiently at least 
when the Weisenberg number Wi, the product of elon- 
gational rate e and polymer relaxation time A, exceeds a 
critical value = eA > 1/2 In simple shear flow 

the ratio between the symmetric and the antisymmetri- 
cal parts of the velocity gradient tensor is one. The rota- 
tional part leads to a tumbling movement of the polymers 
|l^ and the stresses introduced by the elongational part 
are averaged over all orientations. In pure elongational 
flow the polymer stress can grow until the macroscopic 
elastic stress is orders of magnitude larger than the vis- 
cous part. Consequently, the elongational viscosities of 
our PEG solutions measured with a capillary break up 
extensional rheometer (CABER, Thermo Electron Ger- 
many, Karlsuhe) is in the order of several Pas, compared 
to T^e = 3 X 7.2mPas for the solvent. Frequency depen- 
dent measurements of the complex elongational viscosity 
r]l do not exist for our solutions, but for vanishing small 
amplitudes near threshold it should not differ from the 
complex shear viscosity 7^* in terms of a linear response 
theory. 

We now can understand why the polymers do not affect 
the critical accelerations: It is a particularity of surface 
(or interfacial) waves that beside a small viscous sublayer 
near the interface the flow profile can be described by a 
potential and the flow is rotational free and purely elon- 
gational. The elastic part of the linear complex viscosity 
T]'' of our solutions is too small and the dissipative part 
T]' does not differ from the solvent viscosity. Only the 
contributions due to the pure elongational character of 
the flow might affect the polymers, but this occurs only 
for Wi > 1/2. 

To estimate when the flow is strong enough to affect 
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FIG. 5: The integrated intensity F/ of the temporal and an- 
gular averaged spatial Fourier spectrum (Fig. 2J of the inter- 
facial states versus the driving acceleration. The lower liquid 
is the water sugar solution (circles), a 50 ppm PEO solution 
(squares) or the 110 ppm PEO solution (triangles). 



the polymers we look in a first approximation at only 
one spatial dimension. The velocity potential for a given 
wave number k has the form 



^ = uo /kAsin{kx)cosh{kz)cos{(jot) 



(2) 



with A the deformation amplitude of the wave and z 
the direction of gravity. At the interface z = the flow 
is strongest and we can set cosh{kz) = 1. It follows 
for the amplitude Al oi the velocity gradient tensor and 
thus the maximum elongational rate that occurs at the 
position of maximum peak elevation and at the inter- 
face Cmax = Al = kcoA. Only at a driving strength of 
ao = 40m/5^ the amplitude of the basic wave number 
is A = 1mm and the elongational rate is e = 200 1/s 
and Wi = 1.5. At this driving strength the amplitude 
of the basic mode is strong enough to affect the poly- 
mers and we see differences in the amplitude of Fa at 
ko, compared to the pure sugar water solution. The in- 
tegrated Intensity Fj is dominated by the power of the 
fundamental mode and the effect can bee seen best in 
Fig. O However, already at low driving strength the 
wave number continuum at higher wave numbers adds to 
the velocity gradient tensor. The polymers get stretched 
by the high frequency fluctuations and their amplitude is 
reduced (Fig. Apparently, these fluctuations are the 
reason for the disorder transition. 

A comparison with drag reduction in fully developed 
turbulence is instructive. According to Landahl Q the 
occurence of turbulent eddies in wall bounded systems 
like pipes is related to bursting motions, the so-called 
streaks, regions of flow with a high elongational com- 
ponent. The streaks occur near the wall and interact 
subtle with the turbulent eddies in which the energy is 



dissipated via a cascade to smaller lengths scales. The 
pronounced elongational viscosity is supposed to inhibit 
the streaks, yielding less eddies and less turbulent dissi- 
pation. Experimental findings support this picture and 
even a direct correlation of elongational viscosity and the 
ability of a polymer to perform to turbulent drag reduc- 
tion has been found 15]. 

In conclusion we have measured and described the ef- 
fect of a small amount of a flexible polymer on the tran- 
sition from an ordered to a weak turbulent state in para- 
metrically driven inter facial waves. Without polymers we 
observe a rapid transition from a square to a disordered 
state with increasing driving strength. The addition of 
the polymers suppresses the high frequency fluctuations 
and stabilizes the ordered pattern. We find that the pro- 
nounced elongational viscosity is responsible for this ef- 
fect and we suppose that further investigations on the 
effect of polymers on weak turbulent flow might help us 
to understand some of the key features of the scientific 
and technical important problem of turbulent drag re- 
duction, too. 
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